How particle collisions increase the rate of accretion 
from the cosmological background onto primordial black holes 

in braneworld cosmology 
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It is shown that, contrary to the widespread opinion, particle collisions considerably increase 
accretion rate from the cosmological background onto 5D primordial black holes formed during the 
high-energy phase of the Randall-Sundrum Type II braneworld scenario. Increase of accretion rate 
leads to much tighter constraints on initial primordial black hole mass fraction imposed by the 
critical density limit and measurements of high-energy diffuse photon background and antiproton 
excess. 
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It is well known that primordial black holes (PBHs) 
open up a way to obtain information on earliest era of 
the Universe [l |2| . Recent works 0, S H H 0, show 
that PBHs can also be a probe of braneworld cosmolo- 
gies. The Hawking temperature drop and consider- 
able PBH mass growth due to radiation dominated back- 
ground accretion occur 4] in the simplest braneworld cos- 
mology of the Randall-Sundrum Type II (RS2). These 
modifications of PBH properties provide new constraints 
on the initial PBH mass fraction 0, 0- PBH mass 
growth is exponentially sensitive to the accretion effi- 
ciency. However, authors of 0, considered the 
accretion efficiency as a free parameter. In fact, they 
started with F equals 1 in the case of collisionless radia- 
tive background and assumed that it can only decrease 
when relativistic particle scattering length is small. They 
also predicted that a sufficiently low value of F leads to 
weaker constraints than those obtained in the standard 
4D cosmology. It also allow the restriction of the RS2 
curvature radius I from below. Meanwhile, according to 
the theory of spherical accretion of nonrelativistic parti- 
cles 0, the hydrodynamic accretion theory should 
be used instead of collisionless accretion theory [ill IT^ 
for small scattering length. Since "the presence of col- 
lisions between particles restricts tangential motion and 
funnels particles effectively in the radial direction for effi- 
cient capture " 3 , the hydrodynamic accretion rate sub- 
stantially exceeds the collisionless one in the case of non- 
relativistic particles. The aim of the present paper is to 
demonstrate that collisions of relativistic particles anal- 
ogously increase the accretion rate of the radiation dom- 
inated background and lead to the new, much tighter, 
constraints. 

Let us remind the reader that in the RS2 model our 
Universe is a 4D Lorentz metric hypersurface called a 
brane, embedded in a Z2 symmetric 5D anti-de Sitter 
spacetime called the hulk, characterized by the curvature 



radius I, presently bound by condition I < 0.2mm. 
A black hole (BH) behaves as an essentially 5D object 
described by the metric 



dr^ 



r^dnl, 



where dft2 is the line element on a unit 2-sphcrc, if its 
modified Schwarzchild radius 




(2) 



is much smaller than I. We use the units where c = 
G — 1 and introduce 4D Planck scale quantities I4, M4 
and ^4 for convenience. The high energy brane regime is 
characterized by scale factor, energy density, horizon 
mass, and duration 
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where t is cosmic time, and t^ < tc. It is usually assumed 
[1I1I1I10000 that at formation time t, the PBH 
mass is 

= fMh{t,). (4) 

The condition rsH <C I is fulfilled if / ~ 1 and ti < tc- 
The accretion rate 



dM 
~dt 



= Fiirr 



BH 
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and energy density (3) lead to exponential BH mass 
growth 
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at i < 0- In fact, Eq. (5) is justified only in the 
case of collisionless relativistic particles when accretion 
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efficiency F equals 1. However, a simple comparison of 
the scattering length Igc with the radius (2) demonstrates 
that the cosmological background cannot be considered 
as a coUisionless particle medium in the case of PBH 
masses most relevant for observational constraints. In 
fact, the opposite condition Isc ^ rsH is true, in gen- 
eral, for such PBHs. It allows one to treat the accret- 
ing cosmological background as a continuous medium. 
Therefore, to obtain F one should evaluate the accretion 
rate (5) in hydrodynamic approximation. 

Since radiation dominated background is characterized 
by an indefinite particle number, the momentum conser- 
vation equation 



(7) 



should be used together with the entropy conserva- 
tion equation 



('Tw'=);fc = 0. 
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Here comma and semicolon denote usual and covariant 
derivatives; p, ct, P and are energy and entropy den- 
sities, pressure and 4-velocity of accretion flow, respec- 
tively. In the metric (1) adiabatic accretion equations (7) 
and (8) take the form 



Solving Eq. (13) and (14) for u' and T', one obtains 
(compare with 9]) 



where 



D 



and 
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To avoid singularities in the accretion flow, the solution 
of Eqs. (15) must pass through a "critical" or "sound" 
point r = Vg where D = Di = D2 = 0. Solving equations 
D = Dl = D2 = 0, one obtains 



V2 



(19) 



and 



uu + 



' BH 



and 
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BH 
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in light speed units. 

Integration of Eq. (13) leads to a conservation law 
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where u = is the radial component of 4-velocity and 
primes denote the derivative with respect to r. The ra- 
diation dominated background is described by equations 

m 



p{T) = 3P(r) = ^9effT^ 



and 
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where gef / is the number of relativistic particle species 
at temperature T. The value 5e// — 100 will be used 
below because T > ITeV a.t t < tc- Equations (11) 
and (12) allow to rewrite Eqs. (9) and (10) in terms of 
temperature 
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where is the background temperature defined by ex- 
pressions (3) and (11) for energy density 0. Substituting 
Eqs. (19) and (20) into Eq. (21) one obtains 



(22) 



Direct integration of Eq. (10) shows that the term 
r'^a{r)u{r) does not depend on r, while Eqs. (19), (20) 
and (22) allow to evaluate its value 

r'^<j[r)u{r) — r1(j{rs)u{rs)- (23) 

Applying the latter equation sufficiently far from BH, 
where T(r) n, and using Eqs. (12), (19), (20), and 
(22), one obtains 
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' BH 
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r > rsH- 



(24) 



Substituting r > rsn, T{r) ~ Tf,, p{r) ~ p{Tb) = p{t) 
and Eq. (24) into the formula 
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we finally come to the accretion rate (5) with the accre- 
tion efficiency 



F ■ 



3\/3 



1.54 



describing radiation dominated background accretion in 
hydrodynamic approximation applicable in the limit of 
frequent particle collisions. Since (26) exceeds 1, particle 
collisions, in fact, increase both accretion rate (5) and 
PBH mass growth coefficient (6). 

Let us now demonstrate that continuous approxima- 
tion is indeed valid for the most of the high-energy 
phase. Since the characteristic space scale (19) of ac- 
cretion flow is determined by tbhi one should compare 
tbh with the scattering length = lsc{t) ~ l/an. Here 
n = n{Tb) ~ p{Tb)/3Tb is the particle number density 
and (7 = cr(Tfc) ~ cP' /T'^ is the scattering cross section in 
which a = a{Tb) is the interaction constant. 

Let us remind the reader that a substantial mass 
growth will be experienced by PBHs with initial masses 



The lower limit 
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follows (3|] from the absence of excessive large-angle mi- 
crowave anisotropy while the upper one 



max ~ ^ 'j~ ' ^'^'^ 
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follows from both the condition tbh = ^ and Eq. (4) at 
t = tc- Note that BHs with larger masses and radii are 
the usual 4D BHs which do not experience considerable 
accretion mass growth 0, Q . 

For simplicity we will completely neglect scattering 
taking _F = 1 if rsH < he and treat the radiation 
dominated background as a continuous medium taking 
F ~ 1.54 if vbh > he- Let us introduce the ini- 
tial PBH masses Mn and Mi2 for which the condition 
rBH = Isc fulfills, respectively, at the final moment t — 
of the high-energy phase and at the moment t = ti oi 
the PBH formation. These masses allow us to separate 
the mass intervals corresponding to different scenarios 
of PBH mass growth (see Fig. 1), namely, the interval 
Mmin < Mi < Mil where Isc > i"bh and F = 1, the in- 
terval Mi2 < Mi < Mjnax where Isc < rBH and F ~ 1.54 
and the mixed interval Mn < Mi < Mi2 for which a mo- 
ment ti < ti, < tc exists at which rBH{t*) — lsc{t*) and 
one, correspondingly, has tbh < he, F — 1 a.t ti < t < t^ 
and Tbh > he, F ~ 1.54 at t* < t < tc- Equation (6) can 
now be directly applied to each interval with constant F 



giving a composite formula 
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Mil < Mi < M,2; 

ax 

(30) 

for the coefficient of accretion mass growth during the 
high-energy phase. The solid lines on Fig. 1 show the co- 
efficient of accretion mass growth as a function of Mi for 
curvature radii l/h = lO^^lO^^, and lO^^ and / = 0.1. 
Long dashed and short dashed lines correspond, respec- 
tively, to the values F = 1 and F = 0.5 used in 0,13, 
From Fig. 1 one can see that Eq. (30) predicts consid- 
erably larger mass growth than Eq. (6) with F = 1 and 
0.5. One can also see that the region Mj2 < Mi < Mmax 
in which Eq. (30) reduces to Eq. (6) with (26) is quite 
wide. 
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FIG. 1: Accretion mass growth coefficient vs the initial PBH 
mass for choices l/li = 10^^,10^^ and 10^^ (from bottom to 
top). Solid lines correspond to the PBH mass growth co- 
efficient (30), while long dashed and short dashed ones are 
obtained from Eq. (6) with fixed accretion efficiencies F = 1 
and F = 0.5, respectively. 



Let us keep in mind that the curvature radius is experi- 
mentally bound from above by the value 0.2mm ^ 10*^^/4 
[isj . while its lower limit is determined by the choice of 
5D BH evaporation time 0,1111111111 
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where g = 0.023. Since the most direct ways of PBH 
search are connected with PBHs completing their evapo- 
ration now or at a future time, we limit I by Imin = 
which follows from (31) at tevap ~ to — l3.7Gyr and 
M{tc) = Mmax- Note that at l/U < 10^°, one wiU not 
see any difference from the 4D case. 

The accretion efficiency (26) allows us to tighten the 
constraints imposed on the initial PBH mass fraction 0, 
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where ppBH,Mi i^i) mass density of PBHs with for- 

mation mass on the order of Mi and PradiU) is the ra- 
diation energy density at the moment ti of such PBH 
formation. The mass Mi should be evaluated substi- 
tuting the mass M(tc) found from Eq. (31) at given 



I and te 



to into Eq. (30). The PBH mass frac- 



tion grows with time starting from the initial value (32) 
proportionally to a{t)M{t) aX ti< t < tc and propor- 
tionally to a{t) at tc < t < to 0. Describing the PBH 
mass growth by Eq. (30) and restricting the present PBH 
mass density, by 30% of the critical density one obtains 
new constraints on presented by the solid line in Fig. 
2. The latter are much tighter than both the 4D con- 
straint ai < f~^/'^10~^^ and the old 5D constraints 6] 
illustrated by long dashed and short dashed lines evalu- 
ated using Eq. (6) with F = 1 and F = 0.5, respectively. 
The dotted lines both in Fig. 2 and in Figs. 3, 4 below 
are evaluated using the third line of Eq. (30) in order 
to illustrate the degree of violation of hydrodynamic ap- 
proximation in the region l/U> lO'^" where the accretion 
actually begins in a coUisionless regime and only later 
switches over into a hydrodynamical one being, thus, ad- 
equately described by the second line of Eq. (30). 
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FIG. 2: Constraints on initial PBH mass fraction imposed 
by the total matter density. Dotted, long dashed and short 
dashed lines are obtained using PBH mass growth coefficient 
(6) with F = 1.54, 1 and 0.5, respectively. 

Equation (30) also allows to revise the constraints fol- 
lowing from the measured high-energy diffuse photon 
spectrum 0| and antiproton excess Q . Fortunately, it is 
sufficient to consider the PBHs intensively evaporating at 
present in both cases 0,I3- The mass M*{1) = M{tc) oc 
l^^/^ and Hawking temperature T^^ = l/2TTrBH oc 
y/M*{l)l cx r^/-* of such PBHs can be found substi- 
tuting tevap = to into Eq. (31) and then M*[l) from Eq. 
(31) into Eq. (2). The Hawking temperature determines 
the photon energy E ^ bTg^^ with b ~ 5 W\ increasing 
proportionally to l^^^'^ from 250 keV to 150 MeV with 



the curvature radius variation from I/I4 = 10^^ to 10^° 
13 . The diffuse photon spectrum can be estimated 0, Q 
by the formula 



HE) 



M* 



4tt E 

where c is the speed of light and 



n{to,M,{M*)) 



(33) 



n(io,M,(M*)) 



217/4^^9/4 a^{to 
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is the present number density of PBHs with mass of 
the order M*, where teq = 72.6kyr 0, 13 is the time 
of matter-radiation equality and the mass Mi should be 
found by substituting M{tc) = M* into Eq. (30). In par- 
ticular, one can find that initial masses of PBHs complet- 
ing their evaporation at present are equal to Mi = 10"'' 
and lO^g for I/I4 = lO^-^ and 10^°, respectively. Compar- 
ing (33) with the observational data collected in j7|, one 
obtains new constraints presented in Fig. 3 by a solid 
line. These constraints are much tighter than both the 
4D constraint ai < 10~^^ and the old 5D ones, presented 
by dashed lines. 
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FIG. 3: Constraints on initial PBH mass fraction imposed 
by the high-energy diffuse photon background. Dotted, long 
dashed and short dashed lines are obtained using Eq. (6) with 
F = 1.54, 1 and 0.5, respectively. 

The antiproton fiux from brane PBHs can be schemat- 
ically written in the form 



dM, 



oc 



dM* M,M* 



Gn{to,M,{M*)), 



(35) 



where MceV = Mcevil) is the mass of a PBH with 
Hawking temperature Tbh ^ IGeV and the factor 
G ~ 10^ represents the galactic matter density enhance- 
ment at the solar neighborhood. Equation (35) differs 
from equations of ref. Q , combined together in a simi- 
lar way, solely by the relation (30) of M* = M{tc) with 
Mi. This "small" difference, however, allows us to con- 
siderably restrict the " most important result" of Ref. [3 
concerning the p fiux dependence on I. Indeed, using the 
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PBH mass growth coefBcient (6), one obtains 
(35) that $p c>c 1^ with exponent 



P - 



40F - IStt 
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from 



(36) 



Since p < at F < 1 it was argued in Q that the p 
flux is a decreasing function of I (see both dashed hues 
in Fig. 4) and should be used to set a lower bound on 
I. However the new value F ~ 1.54 leads to the positive 
exponent p ~ 0.81 corresponding to increasing p flux (35) 
and decreasing constraints on Ga; (solid line in Fig. 4) up 
to I/I4 ~ 10"^°, where the hydrodynamic approximation is 
violated. These constraints are again much tighter than 
both the 4D constraint a; < 10~^^ and the old 5D ones, 
presented by dashed lines. 
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Because of violation of hydrodynamic approximation 
at I/I4 > 10^° and in the corresponding photon energy 
region E < dOOkeV all the considered constraints reach 
their maxima at I/I4 ~ 10^°. Figures 2-4 demonstrate 
that the tightest constraints following from the critical 
density limit and measurements of the high-energy diffuse 
photon background and antiproton excess at Z/Zq ~ 10'^° 
amount, respectively, to ai ~ 10"'^^, 10""^^ and 10"'^^. 
Since the tightness of all the constraints increases with 
the curvature radius I, all of them can be used to set an 
upper bound on I at I/I4 < 10"^°. 



To summarize, particle collisions make hydrodynamic 
approximation adequate to describe the accretion of the 
radiation dominated background onto PBHs during most 
of the high-energy phase of the RS2 braneworld scenario. 
Continuous background is characterized by 54% higher 
accretion efficiency onto PBHs than that of coUisionlcss 
relativistic particles. This higher accretion efficiency pro- 
vides larger PBH mass growth during the high-energy 
phase and makes it possible to obtain much tighter con- 
straints on initial mass fraction than those found in 
0, 0, 111 as well as than the constraint equal to 10^^^ 
obtained in the standard 4D cosmology. The tightness of 
the obtained constraints mostly increases with curvature 
radius allowing to restrict it from above. 



FIG. 4: Constraints on initial PBH mass fraction imposed by 
the antiproton excess. Dotted, long dashed and short dashed 
lines are obtained using Eq. (6) with F = 1.54, 1 and 0.5, 
respectively. 
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